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ABSTRACT: A fatigue damage computational algorithm utilizing a multiaxial, isothermal, con- 
tinuum-based fatigue damage model for unidirectional metal-matrix composites has been imple- 
mented into the commercial finite element code MARC using MARC user subroutines. Damage 
is introduced into the finite element solution through the concept of effective stress that fully 
couples the fatigue damage calculations with the finite element deformation solution. Two ap- 
plications using the fatigue damage algorithm are presented. First, an axisymmetric stress analysis 
of a circumferentially reinforced ring, wherein both the matrix cladding and the composite core 
were assumed to behave elastic-perfectly plastic. Second, a micromechanics analysis of a fiber/ 
matrix unit cell using both the finite element method and the generalized method of cells (GMC). 
Results are presented in the form of S-N curves and damage distribution plots. 

KEYWORDS: continuum damage mechanics, cracking, fatigue (materials), fracture (materials), 
metal-matrix composites, coupled deformation damage, uncoupled deformation damage, finite 
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In advanced engine designs, materials that allow higher operating speeds and longer dura- 
bility in addition to decreased weight are desirable. The use of metal-matrix composites 
(MMCs) may provide these benefits. For example, titanium metal-matrix composite (TMC) 
rotors are projected to have significant benefits in terms of increased rotor speeds and lower 
weight, as compared to the nickel and titanium rotors currently in service. However, to fully 
realize the benefits offered by MMCs, computationally efficient design and life prediction 
methods must be developed. Analysis of typical aerospace structures subjected to complex 
thermomechanical load histories requires the use of computational approaches such as the finite 
element method. In this regard, it is desirable to develop a life prediction algorithm that can 
be used in conjunction with the finite element method. 

Historically, two basic approaches have been used in predicting the life of structures; un- 
coupled or fully coupled deformation-damage methods. A typical uncoupled analysis consists 
of obtaining the stress state for each element from a finite element analysis, and then, using 
the stress state data as input to a fatigue damage model, the number of cycles to the initiation 
of a crack are predicted. Subsequently, a “local” fracture mechanics approach is then used to 
propagate the crack. For example, a new finite element mesh is constructed to model the crack- 
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tip zone using a series of double nodes. The propagation of the crack is then controlled by a 
strain energy release rate criteria in conjunction with a node release scheme. 

An alternative to the uncoupled method is a fully coupled deformation and damage method. 
By utilizing the concept of effective stress/strain equivalence, the effects of damage are ac- 
counted for in the finite element solution. Through the degradation of the material properties, 
as will be described in the computational scheme section, for a specified increment in damage, 
the individual finite element material properties are degraded and the subsequent finite element 
analysis then calculates the corresponding stress redistribution caused by the damage. This 
sequence of material degradation (representing damage) followed by re-analysis to capture 
stress redistribution effects is repeated until structural failure. 

This is the approach taken in the present study. Specifically, the computational scheme 
developed uses MARC, a nonlinear finite element code in which the fatigue damage algorithm 
is coupled to MARC through the use of provided user subroutines. In the next section, the 
requisite fatigue damage equations are presented. In the following section, the computationally 
coupled fatigue damage algorithm will be outlined. Finally, two applications ot the computa- 
tional scheme will be presented. The first example is a reinforced MMC ring, representing a 
typical engine component, and the results will be presented in terms of the evolution of damage 
in the ring cross section. The second example is a micromechanics analysis ot a fiber-matrix 
representative volume element (RVE). The damage distribution in the matrix will be shown. 


Fatigue Damage Formulation 

The fatigue damage calculations utilize a recently developed multiaxial, isothermal, contin- 
uum damage mechanics model for the fatigue of unidirectional metal-matrix composites |/|. 
The model is phenomenological, stress-based, and assumes a single scalar internal damage 
variable. Note that for an initially anisotropic material, the evolution of the damage, although 
a scalar, is directionally dependent. As w ill be shown, this directional dependence is accounted 
for in the terms, F m , <t>, h and <1>„. The present multiaxial, isothermal, continuum damage model 
for unidirectional metal-matrix composites may be expressed as [ / 1 
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N is the number of cycles at the current stress state (a k ) and the increment in damage, (D k ~ 
D k _ |), where D k and D k , is the amount of damage at the current and previous increments, 
respectively. The quantity, a, that is a function of the current stress state is defined as 


a — 1 - a 
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where ( ) are the Macauley brackets. In Eq 2, the fatigue limit surface, <Jy and the static 
fracture surface, <F W , are defined as 
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Note, the case, = 0, indicates static fracture, which is failure, making it unnecessary to 
perform the fatigue calculations. Thus, having to consider the possibility of a being undefined 
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is unnecessary. As will be discussed in the following section, the finite element is considered 
to have failed completely. The case, (<f> /7 ) - 0, indicates that the current stress state is below 
the fatigue limit and thus a is set equal to 1. This presents a special case when integrating the 
fatigue damage expression, Eq 1, and will be considered later in this section. 

The quantity, F„ n the normalized stress amplitude, is defined as 


- 1 max max 
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In the preceding equations, /<, is the time at the beginning of the current load cycle, and t is 
some time during the load cycle. The general form for may be expressed as 
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It is here in Eq 6 where the evolution of the damage becomes directionally dependent. This 
simply amounts to the assumption of partial anisotropy, where the “extent” (magnitude) of 
damage is affected by the directionality of the stress state. Specifically, the directional depen- 
dence enters through the quantities, /, , / 2 , / 3 , w ( ,, and r) { >. The quantities, , I 2 , / 3 , are invariants 
having the form 


7 > = \ - d t djS jk S kl + (djdjSijf 

I 2 = dJjSjj'Su - ( djdjSij ) 2 (7) 

/ 3 = (dfdjSjj) 2 

that are a function of the current deviatoric stress state, Sj, as well as the vector, d h that defines 
the materials fiber orientation. In addition, the terms, , and if t represent the ratios of 

longitudinal to transverse normal and shear stresses, respectively. Note, the longitudinal direc- 
tion is parallel to the fiber direction and transverse is perpendicular to the fiber direction. For 
initially transversely isotropic materials, o> { > and 17, } are > 1, and for isotropic materials, , 
and rf ( , are = 1. 

For a current state of stress, c T k , which is above the fatigue limit, that is, a ^ 1 (integrating 
Eq l ) results in an expression for the number of cycles, N, that is 

(M - (1 - D k f + T '" - [j - (1 - D,-,)^ 1 ] 1 -") 

F p m (\ - a)((3 + I) (8) 

Note that D*_, is the total amount of damage at the beginning of the load block and D k is 
the total amount of damage at the end of this load block. 

In the present computational scheme, since the damage increment is controlled, both D k and 
D k _ 1 are known. That is, D k = D k j T- AD where AD is the user-specified increment in damage. 
Thus, Eq 8 is used to predict the increment in the number of cycles for each element, N e , due 
to the increment in damage. 

To calculate the number of cycles to failure for an initial damage amount, D*_,, let D k = 1, 
which results in the following 


(1 ~ [1 ~ (1 ~ D k - if -'] 1 -") 
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As will be shown in the following section, it is also necessary to rewrite Eq 8 in terms of 
the damage D k , that is 


lh = 


(I - D k l )(' M T 
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Now, consider the case in which the current stress state is below the fatigue limit, that is, 
(<!>„) = 0, which leads to a k 1. Thus, Eq 1 takes the form 
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Upon integrating Eq 11, the increment in cycles, N, with initial damage, D k ,, may be expressed 
as 


N = 


log[ 1 - (1 - D k f" | log| 1 ~ (1 ~ D k ,r"i 
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Alternatively, the following expression for the damage, D k , may be expressed as 
D k = I ~ {1 |1 - (1 - D k ,/ <M | exp((/3 + l)T^)l}/o i 

For the number of cycles to failure, let D k — 1 
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The effect of damage is included in the finite element stress analysis by utilizing the concept 
of “effective stress" [2] and the hypothesis of strain equivalence [J,4J. Thus, the effect of 
damage may be accounted for by either using the effective stress, that is, increasing the stress, 
or by simply degrading the elastic and plastic material properties to represent the material 
softening due to damage. Material softening is used in the present approach with the degraded 
elastic constitutive matrix calculated by 

[C| = (1 - D k )[C] (15) 


and similarly, the plastic material properties, for example, yield stress, a v , is given by 


oy = ( l - D k )cr y 


(16) 


Computational Scheme 

General Framework 

The present version of the fatigue damage algorithm utilizes average quantities in the damage 
calculations. For example, the stresses for each integration point are determined and then all 
integration points are averaged to give one stress state for each element. All subsequent damage 
calculations use these average quantities. However, the program was written in a sufficiently 
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general form so that all of the damage calculations may be performed at each integration point 
with minimal modifications. Specifically, all that is required is increased dimensions for various 
sorage arrays. Note that, because of the scalar nature of the damage parameter used here, 
together with the assumptions previously discussed, the task of local (constitutive) level inte- 
gration is simply reduced to the use of the “exact” closed-form expressions given in the 
preceding section. b 

A flowchart of the developed life prediction scheme is shown in Fig. 1. The deformation 
analysis is the actual finite element run. First, note there are two levels of failure criteria checks- 
element level and structural level. The element level includes a static fracture surface check 
that is part of the fatigue damage model (Eq 4). Additional criteria may be included such as a 
c eck on total mechanical strains, etc. If an element violates one of these failure criteria that 
element is considered to have “failed” and its damage, D, is set equal to the maximum amount 
of damage allowed. As shown in Fig. I, for a coupled analysis, the damage calculations are 
emunated and a deformation analysis is rerun in order to account for the stress redistribution 



-► Uncoupled Analysis 
Coupled Analysis 

FIG. 1 Coupled/uncoupled life prediction scheme. 
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due to that element's failure. For an uncoupled analysis, when element lailure occurs, the 
analysis does not loop back and perform another deformation analysis; instead, it continues to 
the next element and performs the damage calculations. The structural level criteria monitors 
the global response of the structure. This could take the form of a check upon selected noda 
displacements that, if they violate a specified displacement criteria, the structure is considered 
to have failed. For example, the tip displacement of a turbine blade may be required to stay 
within a given tolerance. Again, note that for an uncoupled analysis no structural failure criteria 
check is performed since the present fatigue damage algorithm assumes that the structure, when 
subjected to the initial “applied” load cycle, is in a completely undamaged state. 

The applied load cycles just mentioned are those actually defined by the user through t e 
finite element input data file. In addition, the number of load increments per applied cycle must 
be specified by the user. This is necessary so that the program can internally monitor when a 
given applied load cycle has been completed and begin the damage calculation phase. The 
number of load increments used in the load cycle usually depends on the nonlinearity ol the 
structural response and requires experience on the part of the user. At each increment during 
the applied load cycle, the average stress (strain is optional) state for each element is stored. 
At the last load increment of the current applied load cycle, the DAMAGE subroutine is called 
to perform the fatigue damage calculations, see Fig. 2. 

When DAMAGE is entered for each element, various element quantities, such as a, 4>„, <1>,„ 
p are calculated and stored. When DAMAGE has been called for the last element in the 
mesh, the fatigue damage calculations are performed. Figure 2 shows the general algorithm tor 

the fatigue damage calculations. , 

Presently, the damage calculations are controlled by the increment in damage. \D. I he user 
specifies the allowable increment in damage, for example, AD = 0.15, (15%). In CALCN, 
using Eq 9 or Eq 12, based upon the new value of damage and the given element's stress state, 
the number of cycles to failure, N' h , is calculated and stored. Next a “sorting” subroutine, 
SORTN, is called (see Fig. 2) to determine which element has the minimum number ol cycles 
to failure and is chosen as the “controlling element,” that is 

N hlm „ = min/V, e = 1 -» numel (17 > 



FIG. 2 — Fatigue damage calculations . 
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Once the controlling number of cycles has been determined, the corresponding, actual 
amount of damage, /)*, in all of the remaining elements must be recalculated. This is performed 
in subroutine CALCD using Eq 10 or 13. Note that since the damage was incremented by a 
specified amount, the controlling element’s damage is already known. 

Figure 3 shows the cycle scheme used in the code. Recall that the “applied’ ’ load cycle is 
the actual load history that is applied and used in the finite element analysis. The subsequent 
cycles, shown in dashed lines, are the predicted cycles corresponding to N Fmm that is determined 
in SORTN. Here, it is assumed that the stress state in each element remains constant during 
the predicted N Fmin cycles and at the end of A V mjn , each element has incurred an amount of 
damage as calculated in CALCD. Note for a coupled analysis, the next applied load cycle is 
run in the finite element analysis to account for the stress redistribution due to the new damage 
state in each element (that is, lY k ) and again a new /V / n , in is predicted. This sequence of applied 
load cycle and predicted cycles is repeated until the structure has failed. For an uncoupled 
analysis, only one sequence, that is, only one applied load cycle, is performed. The resulting 
AVmin would be used to merely indicate the location of damage initiation. 

In preparation for the subsequent applied load cycle, the element material properties are 
degraded according to the newly calculated element damage, D‘ k . Finally, a subroutine, 
PATSTR, was also written that generates PATRAN [5] element results files. These files contain 
damage distributions at specified increments during the fatigue damage analysis. In addition, 
output files containing a summary table of the current number of fatigue cycles and remaining 
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cycles to failure for each element, and a summary table showing the damage evolution in each 
element are also generated. 

Some general comments on the fatigue damage algorithm need to be made. First, once an 
element attains the user specified maximum allowable damage, the element is assumed to fail 
and is no longer considered in any subsequent damage calculations. In the example presented 
in the next section, a cutoff value of 95% was specified based upon preliminary experience 
with difficulty in achieving global convergence when the element stiffness was reduced below 
5%. Further investigation of convergence difficulties needs to be addressed. Second, in the 
present fully coupled damage-deformation analysis, a perfect plasticity model was used. This 
idealization eliminates the need to account for cyclic hardening, and the corresponding update 
of the internal variables, which may occur during a specific block of fatigue cycles. In order 
to accurately account for the hardening, a projecli on/update of the internal variables through 
the load block would be required. In addition, it is usually assumed that the fatigue damage 
calculations are applied to a “stabilized” stress redistribution. Thus, when hardening is present, 
more than one applied load cycle may be necessary in order to achieve the stabilized stress 
redistribution. 

MARC Implementation 

As an example, the present fatigue damage algorithm was implemented into the finite element 
code MARC. Please note that based upon the discussion in the previous section, it should be 
apparent that the fatigue damage algorithm was developed in a form that is independent of the 
specific finite element code that is used. Most commercial finite element codes provide capa- 
bilities similar to those described here with regards to MARC. 

MARC provides various user subroutines [6] that allow implementation of constitutive mod- 
els, failure criteria, new elements, etc. By using a few select MARC user subroutines, the 
continuum-based fatigue damage model has been coupled with the nonlinear finite element 
solution scheme. The primary MARC user subroutine required is ELEVAR. In addition, when 
using the plasticity model contained in MARC, the subroutines HOOKLW and ANPLAS are 
used to degrade the elastic and plastic material properties, respectively. In this study, the user 
subroutine HYPELA was also used to implement the Bodner-Partom viscoplastic model into 
MARC. 

The subroutine ELEVAR is called at the end of each load increment once global convergence 
has been attained, and is intended to be used to output element quantities at the end of a given 
increment. In this algorithm, ELEVAR is used to store the current converged stress state for 
each element during the applied load cycle. The meaning of applied load cycle is the same as 
discussed earlier in the previous section. This is done through the MARC subroutines 
HOOKLW, for the elastic constants, and ANPLAS, for the anisotropic yield stress ratios using 
Eqs 17 and 18. Again, HOOKLW and ANPLAS are used when the MARC plasticity model is 
used. For the Bodner-Partom viscoplastic model, similar calculations are performed within the 
HYPELA routine. 

Example Applications 

A Cladded MMC Ring Insert 

As stated previously, one of the primary motivations of this research is to establish a com- 
putationally efficient method for predicting the fatigue life of typical aerospace components. 
This includes the ability to predict the location of damage initiation and to be able to track the 
propagation of damage throughout the structure. With this in mind, the fatigue damage algo- 
rithm was applied to a cladded MMC ring. The reasons for choosing this specific structure are 
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TABLE 1 — Material properties for plastic deformation model 
(MPa). 


SiC/Ti 15-3 composite material 

Elastic 

= 183959, E 2 = E, = 1 14457, v i2 = 0.28, p M - is* - 0.32 
Inelastic 


<r x ' rr v <t, 

MATRIX MATE-RIAL (Ti 15-3) 

Elastic 

E = 74466, v = 0.32 
Inelastic 
rr, = 514 


Noi l — The subscript, 1, denotes fiber direction. 


twofold. First, it represents a MMC rotor insert currently under consideration in advanced 
engine designs. Secondly, because of its axisymmetric geometry and load conditions, qualitative 
stress distributions are known a priori. For example, maximum circumferential stress in the 
core occurs at its inner diameter and likewise for the cladding, thus providing some intuitive 
feel for where damage initiation will occur as well as how it may propagate. 

The composite core was described by Hill's anisotropic elastic-plastic constitutive model 
available in MARC [6|, while the matrix cladding was assumed to be isotropic and elastic 
perfectly-plastic. The elastic and inelastic material parameters required for the deformation 
analysis are given in Table 1, while the associated material parameters for the fatigue damage 
model are given in Table 2. Note that the matrix cladding utilizes the isotropic form of the 
fatigue damage model, that is, to u = (o f7 = (o„, = r) u = -q fl = rj m = 1 , whereas the composite 
core is represented by the transversely isotropic form of the model. The finite element model, 
representing the cross section of the ring, Fig. 4, consisted of 225 nodes and 64 eight-node 
axisymmetric elements (MARC element No. 28). A uniform pressure load was applied along 
the inner diameter of the ring. 

With regards to the deformation analysis, burst pressure predictions have been previously 
made and compared to limited experimental data [7,#]. Very good correlation was observed, 
thus providing a level of confidence in the finite element modeling of the ring. 


TABLE 2 — Material properties for fatigue model |7j (MPa). 


= 10694 
= 1972 
/ 3 = 1.842 
a = 0.012 
M = 22371 


SiC/Ti 15-3 composite material 

a>„ = 5.5 
(o n = 12.482 
w r „ =11.8 
Vu = T//7 = i?,„ 


1.0 


MATRIX MATERIAL (Ti 
(T u = 6081 

(Tfj = 965 
/ 3 = 2.27 
a = 0.0365 
M = 6205 


1 5-3) — isotropic simplification 


= 1.0 
ix)fj = 1.0 

a>„ = L0 


n u = r) fl = 7) m = 1 .0 
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Center 

Line 



fig. 4 — Cladded MMC ring geometry and finite element model: (top) actual ring geometry (di- 
mensions are in mm) and (bottom) idealized ring geometry. 


Two types of fatigue life analyses were performed, namely, an uncoupled and a coupled 
analysis. The uncoupled life prediction results were obtained by taking four elements, labeled 
12, 3 and 4 in Fig. 4 (bottom), in the radial direction of the composite core of the ring. It was 
assumed that the stress state was relatively constant in the ^-direction, thus one element would 
represent all of the elements in a column of the composite core. In the uncoupled analysis, no 
fatigue calculations were performed on the elements associated with the matrix cladding be- 
cause of the initially low stress levels in the matrix cladding causing infinite fatigue lives to be 
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calculated. A finite element analysis was run in which the ring was subjected to a cyclic pressure 
load. During the analysis, the stress state history for each of the four elements was stored The 
stress history for a given element (I, 2, 3, or 4) was then used to predict the corresponding 
taligue life for that element. Since only a single cycle was used, no stress redistribution effects 
due to damage, were accounted for in the uncoupled analysis. Figure 5 (top) shows the results 
ot the uncoupled fatigue damage analysis. As expected, since FJement I has the highest stress 




FIG. 5 — Life prediction for uncoupled (top) and coupled (bottom) analysis. 
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N/N F - 0.20 



FIG. 6 — Fatigue damage distribution m 
1.0. Cycles to failure. N, = 155. 



ring cross section: (left) N/N, = 0.2 and (right) N/N,, ~ 


I, correspondingly has ,he shor.est fariguc life. Thus, damage is predreled to inidare along 

^Tly' analysis follows ,he alg.rifhm presen, ed in Ihe prev.ous 
1 t L^c intn iccount the effects of stress redistribution due to the propagation of 

finite lives and damage are predicted for the cladding 
due to stress redistribution effects. Figure 5 ( bottom ) shows the tully coupled deformation and 
fatigue damage analysis results. The solid circles are analysis results obtained using the present 
fu|iy coupled^ algorithm. A smooth solid line is drawn through these points to obtain the S-N 
[urvV Not the dashed part of the S-N curve indicates an “extrapolation- since no analysis 
w,s run at a PIP less than 0.30. Upon comparison, one observes that at pressures close to 
burst pressure the fatigue life predicted by the coupled analysis is close to that of the uncoupled 
aSs^s "nc; at high stress levels, once the damage initiates in the core, “structural fa, 1- 
ure of the ring occurred rapidly. On the other hand, at low stress levels, the fatigue li e as 
predicted by the coupled analysis is longer than that predicted from the uncoupled analyse 
at the inside^ diameter. This difference may be viewed as the effect of propagation of the damag 
"cross section. This propagation is caused by the stress redistribution effect that 
are automatically captured by performing a fully coupled deformation and fatigue damage 

""Finally Fig. 6 shows two selected damage distribution plots in the ring cross section pro- 
duc^by the coupled fatigue damage analysis for PIP, = 0.85. Note that in Fig^ 6 I left), the 
damage initiates along the inner diameter of the composite core. In 

to completely fail are located in the inner radius “comers of the TMC core. Iheretore 
the assumption used in the uncoupled analysis of a uniform stress state in the z-direction 
for a given column of core elements was not appropriate. Converse y, in tg. 
structural failure of the ring is depicted (that is, the composite core has completely faded 
and due to stress redistribufion, the matrix cladding has accumulated significant amounts of 

damage. 
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FIG. 7 — Finite element mesh (top) and GMC (bottom) RVEs. 


Micromechanical Failure Analysis 

The second application is on the microscale in which a fiber/matrix unit cell representing a 
square pack, 35% fiber volume ratio (FVR), SiC/Ti-15-3 composite was analyzed Two ap- 
proaches were used to model and analyze the unit cell. The first approach utilized a finite 
element representation in which the unit cell was idealized using two mesh discretizations, 
consisting of 4 and 64 eight-node three-dimensional elements. Fig. 7. For the deformation 
response, the Bodncr-Partom viscoplastic model was characterized for the Ti- 15-3 matrix ma- 
terial (see Table 3) and the fiber utilized a simple isotropic linear elastic model (see Table 4). 


TABLE 3 — Material properties for viscoplastic deformation model 
(MPa) SiC fiber material. 


Elastic 

E = 399910, v = 0.25 
Inelastic 

D {) ~ 0, 2o = Z t = m ~ n = Q = 1 





78 APPLICATIONS OF CONTINUUM DAMAGE MECHANICS 


TABLE 4 — Material properties for viscoplastic deformation model 
(MPa) Ti 15-3 matrix material. 


Elastic 

E = 74466, v = 0.32 

Inelastic 

D 0 = 1000, Zo = Z. = 120 
m = 10, n = 3, Q = 1 


The Bodner-Partom viscoplastic model was implemented into MARC using the HYPELA user 
subroutine. The fatigue damage model, taken in its isotropic form, was characterized for the 

Ti- 15-3 matrix constituent. , , . , 

The second approach utilized the generalized method of cells, CMC, developed by Aboudi 
191 GMC is a continuum-based micromechanics model that provides closed-torm expressions 
for the macro response of the composite in terms of the individual constituents (phases) The 
GMC model and the fatigue damage algorithm have been incorporated into the stand-alone 
micromechanics analysis code, MAC 1 10). MAC has the ability to analyze a material volume 
element subjected to various thermal, mechanical (stress or strain control), and thermomechan- 
ical load histories, and wherein different integration algorithms and a variety of constitutive 
models can be selected. Previous work [//] has shown the accuracy and computational elh- 
ciencv of GMC with respect to deformation analysis. Part of this work, and future research, is 
to determine if GMC also provides accurate and computationally efficient fatigue damage 

Pr ' The analysis consisted of stress-controlled cyclic loads applied in the longitudinal direction 
with respect to the unit cell. All four unit cells (both finite clement and GMC) produced identical 



FIG. 8 —Longitudinal stress versus strain response: finite element and GMC. 
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TABLE 5— Fatigue lives for a 35% fiber volume fraction SCS-63Ti-l5-3 composite loaded 

longitudinally. 

Fatigue Life, Cycles 

Finite Element GMC 


Load, MPa 


64 


Difference, % 


4-Cell 


49-Cell 


Difference, % 


1345 

2939 

1207 

4957 

931 

21497 

793 

64496 


2858 

3.0 

4859 

2.0 

20992 

2.4 

62366 

3.4 


2497 2494 

4245 4234 

20459 20098 

61235 59843 


0.2 

0.25 

2 . 

2.3 


longitudinal macro stress-strain responses, as shown in Fig. 8. However, the predicted longi- 
tudinal fatigue lives are different, see Table 5. An example macrolongitudinal S-N curve pro- 
duced with the 64 finite element model is shown in Fig. 9. Note, that if the SCS-6 fiber is taken 
to have a constant ultimate tensile strength (UTS) value, the lower stress amplitude lives in 
Fig. 9 are run outs. Consequently, in order to agree with the experimentally determined fatigue 
lives, a fiber UTS model that degrades as a function of time and temperature is required. Fig. 
10. All the lives indicated in Table 5 were generated using this time-dependent liber UTS 
model. Referring to Table 5, note that the 4 and 64 finite element unit cells give results that 
are consistently within 2 to 3% of each other for all stress levels, and the 4 and 49 cell GMC 
results are w ithin 0.2 to 2% of each other. This shows that for the longitudinal response the 4 
finite element unit cell produces results that are comparable to the 64 finite element unit cell. 
Similarly, the 4-cell GMC representative volume element (RVE) gives results that are as ac- 
curate as the 49-ccll GMC RVE. Table 6 shows a comparison between the 4 finite clement unit 
cell versus the 4-cell GMC RVE and the 64 finite element unit cell versus the 49-cell GMC 



FIG. 9 — Life prediction for SiC/Ti-I5-3 composite. 
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PIG ]0 — Fiber constituent ultimate tensile strength {UTS) failure diagram. 


RVE. Note that at the higher stress amplitudes (1345 and 1207 MPa), the relative differences 
are on the order of 14 to 17%, while at the lower stress amplitudes (931 and 793 MPa), there 
is a relative difference of 4 to 5%. 

Figure 1 1 shows the damage distribution for the fatigue analysis having a maximum stress 
amplitude of 1345 MPa. The damage distributions were taken just before failure of the unit 
cell (that is, fiber fracture). Figure 1 1 shows the corresponding damage distributions tor both 
the 4 and 64 finite element unit cells and the GMC 4-cell and 49-cell RVE models. 


Conclusions 

A coupled/uncoupled deformation and fatigue damage algorithm has been presented. The 
algorithm utilizes a multiaxial, isothermal, stress-based, transversely isotropic continuum fa- 
tigue damage model in which the fatigue damage calculations are coupled with the nonlinear 
finite element solution using the concept of effective stress. Incorporated in the life prediction 
scheme are failure criteria checks at both the element and structural level. The algorithm has 
been applied to a cladded MMC ring insert representing a typical aerospace component and 
results have been presented in terms of S-N curves along the damage distribution plots over 
the ring cross section. The fatigue damage results presented are qualitative in nature since no 


TABLE 6 —Fatigue lives for a 35% fiber fraction SCS-6/Ti- 1 5-3 composite loaded Itmgitudim illy. 


Load, MPa 



Fatigue Life, Cycles 



4 Finite 
Element 

4-Cell 
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Difference, % 
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1345 
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FIG. I \— Representative damage states: finite element and CMC unit cells. 


wTJTnTrf ' reSU I US ^ CU , rren ' ly ava,lahlc However, full-scale burst pressure and fatigue tests 
p rtormcd on similar cladded MMC rings. Once these test results become available a 
similar hnite element analysis will be conducted to verify the present fatigue damage algorithm 
and continuum fatigue damage model. s aigontnm 

With regards to the micromechanics fatigue damage analysis, future work will involve nre 
ctmg the transverse fatigue behavior of the same SiCYTi-15-3 composite. Comparisons will 

of accurT* 1 , etWCen the ,lnite elernent method and the generalized method of cells in terms 
h Tump computational efficiency. Finally, an attempt will be made to repeat the macro 
b d MMC ring insert life ami,*, but this time „ s i„g , m,cromecbu„,c s -b as cCp-oach 
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